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Abstract. 

Water's behavior differs from that of normal fluids, having more than sixty 
anomalies. Simulations and theories propose that many of these anomalies result 
from the coexistence of two liquid phases with different densities. Experiments in bulk 
water confirm the existence of two local arrangements of water molecules with different 
densities, but, because of inevitable freezing at low temperature T, cannot ascertain 
whether the two arrangements separate into two phases. To avoid the freezing, new 
experiments measure the dynamics of water at low T on the surface of proteins, finding 
a crossover from a non-Arrhenius regime at high T to a regime that is approximately 
Arrhenius at low T. Motivated by these experiments, Kumar et al. [Phys. Rev. 
Lett. 100, 105701 (2008)] investigated, by Monte Carlo simulations and mean field 
calculations on a cell model for water in two dimensions (2D), the relation of the 
dynamic crossover with the coexistence of two liquid phases. They show that the 
crossover in the orientational correlation time r is a consequence of the rearrangement 
of the hydrogen bonds at low T, and predict that: 

(i) the dynamic crossover is isochronic, i.e. the value of the crossover time tl is 
approximately independent of pressure P; 

(ii) the Arrhenius activation energy Ea(P) of the low-T regime decreases upon 
increasing P; 

(iii) the temperature T*(P) at which r reaches a fixed macroscopic time r* > tl 
decreases upon increasing P; in particular, this is true also for the crossover 
temperature Tl(P) at which r = tl. 

Here, we compare these predictions with recent quasielastic neutron scattering 
(QENS) experiments performed by X.-Q. Chu et al. on hydrated proteins at different 
values of P. We find that the experiments are consistent with these three predictions. 
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1. Introduction 

Water has many anomalies compared to normal liquids [lj. Experiments show that 
its thermodynamics fluctuations and response functions, such as the isobaric specific 
heat Cp or the magnitude of the isobaric thermal expansion coefficient ap, largely 
increase when temperature is decreased [2]. These anomalies are more pronounced in 
the supercooled liquid state, below 0°C. This state, metastable with respect to ice, can 
be observed at temperatures as low as — 47°C in plants [3J, — 41°C in laboratory at 
atmospheric pressure [I] and — 92°C at 2 kbar [5]. 

The anomalies have been interpreted on the basis of models that propose different 
scenarios. The scenarios can be divided in two main categories: (a) those scenarios 
that include the coexistence at low T of two liquids with different densities, and (6) a 
scenario in which water forms local regions of different densities, but does not separate 
into two phases. 

1.1. Scenarios with coexistence of two liquids 

The liquid-liquid critical point scenario. In 1992 Poole et al. [6], on the basis of 
molecular dynamic simulations for ST2- and TIP4P-water, proposed that supercooled 
water separates into two liquid phases with different densities below 200 K and above 
150 MPa. The low-density liquid (LDL) form of water and the high-density liquid (HDL) 
coexist along a first-order phase separation line with negative slope in the pressure- 
temperature P-T phase diagram and terminate in a liquid-liquid critical point (LLCP). 
The LLCP scenario has been confirmed by simulations with various water models [7j. 

The occurrence of the LLCP has been rationalized with different theoretical models 
by assuming (a) the anticorrelation between energy and volume in the H (hydrogen) 
bonds formation and (b) the possibility of forming different kinds of H bonds. An 
example is the thermodynamic free energy model of Poole et al. that hypothesizes the 
formation of strong H bonds when a geometrical condition on the molar volume is 
satisfied [8]. 

Another example is the microscopic cell model introduced in Ref . [9j \TU\ and studied 
in detail in Ref. [HI [12] , in which the correlation between the H bonds, due to the O- 
0-0 interaction, is incorporated. In this model the liquid-liquid (LL) phase transition 
is due to the tetrahedral ordering at low T and P of the H bonds, as explained in the 
following section. 

The LLCP is predicted to lie in a region of the P-T phase diagram where the 
freezing of bulk water is inevitable. Therefore, direct experimental verification of the 
LLCP scenario is challenging. However, Soper and Ricci in 2000 showed with neutron 
diffraction measurements that the local arrangement of water molecules changes up 
to the second shell, increasing the local density when P is increased from 26 MPa to 
400 MPa at -5.15°C [IB]. This structural change was initially observed by varying T 
from about 263 K to 313 K in x-ray structure factor experiments for heavy water D 2 
in 1983 [14J. The data show that by decreasing T the average 0-0-0 angle increases 
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toward the tetrahedral angle 109.47° [14] . This result has been reaffirmed by Ricci et 
al. in a recent experiment [15j . 

The critical point free scenario. Another scenario that hypothesizes a liquid-liquid (LL) 
phase transition has been considered recently [16]. In this scenario the LL first order 
phase transition extends to negative P and merges the liquid spinodal, but without a 
critical point. A rationalization of this scenario has been recently proposed on the basis 
of the microscopic cell model for water [9[ [17], showing that this scenario has a liquid 
spinodal that reenters from negative to positive P [T7] . 

The stability limit scenario. The reentrant spinodal discussed in Ref. [UJ was earlier 
proposed in the stability limit scenario [18] as the origin of the anomalies of water. 
Although in the initial formulation of the stability limit scenario the occurrence of the LL 
phase transition was not hypothesized, it was successively introduced for thermodynamic 
consistency [19] and found also in the thermodynamic free energy model of Poole et al. 



1.2. Scenario without coexistence of two liquids 

The singularity free scenario. This scenario assumes that the H-bonds linking molecules 
are uncorrelated. Under this hypothesis, water anomalies are the effect of the negative 
volume-entropy cross fluctuations [20j [21] and the large increase of response functions 
seen in the experiments represents only an apparent singularity, due to local density 
fluctuations. The regions with different local density do not form separate phases. 
A pressure increase gives rise to a sharp, but continuous, increase of density, as in 
the supercritical region of the LLCP scenario. The continuous structural change is 
found also in ab initio water simulations [22] at very high pressure, P = 10 4 MPa, and 
T = 600 K. This scenario is recovered by the microscopic cell model [9] in the limiting 
case of no cooperativity among the H bonds. 

The singularity-free scenario, in the region accessible by experiments, reproduces 
the same phase diagram as the scenarios with coexistence of two liquids. Therefore, it is 
interesting to understand if there are oberservable differences among these scenarios. 
In particular, recent experiments and simulations, described in the next sections, 
have analyzed the case of the dynamics of water surrounding proteins or confined in 
nanopores, interpreting the experimental results within the context of the different 
scenarios. 

2. Hydrated proteins 

Recent experiments on surface water and water hydrating proteins [231 EH EE] have 
shown that liquid water exists at temperatures as low as — 113°C [26] at ambient 
pressure. At these extremely low temperatures interesting dynamical phenomena 
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occur [27] EHl ESI ED], suggesting a possible relation for the dynamics of the biological 
macromolecules with that of the surrounding water [25], [31] . 

At low T, proteins exist in a ("glassy") state with no conformational flexibility 
and with very low biological activity. For hydrated proteins above about 220 K, the 
flexibility is restored, allowing the sampling of more conformational sub-states. Hence, 
the protein becomes biologically active at these temperatures. This dynamical transition 
is common to many biopolymers and is believed to be triggered by the strong coupling 
with the mobility of the hydration water [28] , which shows a similar dynamical transition 
at about the same T [23]. Chen et al. by studying the translational correlation time of 
water molecules hydrating a lysozyme protein [23], DNA [32] and RNA [33], found that 
at about 220 K the dynamics of hydration water changes from non-Arrhenius at high 
T to Arrhenius at low T. By definition a correlation time r has an Arrhenius behavior 
when 

r = r exp[E A /k B T] (1) 

where tq is the correlation time in the high-T limit, E\ is a T-independent activation 
energy and hp is the Boltzmann constant. On the other hand, r is non-Arrhenius when 
its behavior cannot be fitted with the expression in Eq. ([1]). 

Motivated by these experiments, Kumar et al. [31] simulated using the TIP5P model 
the dynamics and thermodynamic behavior of hydration water for (i) an orthorhombic 
form of hen egg-white lysozyme and (ii) a Dickerson dodecamer DNA at constant 
pressure P = 1 atm, several constant temperatures T, and constant number of water 
molecules N. Kumar et al. [34] found that the mean square fluctuations (x 2 ) of the 
biomolecules changes its functional form below T p pa 245 K, for both lysozyme and 
DNA. They also found that the specific heat Cp of the total system (biopolymer and 
water) displays a maximum at T\y ~ (250 ± 10) K for both biomolecules. 

To describe the quantitative changes in structure of hydration water, Kumar et 
al. [34J calculated the local tetrahedral order parameter Q [35] for hydration water 
surrounding lysozyme and DNA and found that the rate of increase of Q has a maximum 
at Tq = (245 ± 10) K, the same temperature of the crossover in the behavior of 
mean square fluctuations. Finally, they found that the diffusivity of hydration water 
exhibits a dynamic crossover from non-Arrhenius to Arrhenius behavior at the crossover 
temperature T x « (245 ± 10) K for lysozyme and T x w (250 ± 10) K for DNA. Note 
that T x is much higher than the glass transition temperature, estimated for TIP5P as 
T g = 215K [7]. Thus this crossover is not likely to be related to the glass transition 
in water. Therefore, the fact that T p pa T x pa Tw ~ Tq is evidence of the correlation 
between the changes in protein fluctuations and the hydration water thermodynamics 
and structure. Before analyzing in more details this relation, is worth considering the 
results of experiments about water confined in nanostructures. 
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3. Confined water 

The non-Arrhenius to Arrhenius dynamic crossover has been found also in water confined 
in 20A MCM-41 silica pores. In 2004, indeed, Faraone et al. found the crossover 
for confined water at Tl ~ 221 K by studying the structural relaxation time by 
neutron scattering [36J. This result was reinforced by the neutron magnetic resonance 
measurements of Mallamace et al. showing the crossover at Tl ~ 225 K for the self 
diffusion of confined water |26j. 

In 2005, Liu et al. [37] showed that, by increasing the pressure, the dynamic 
crossover of water confined in 20A MCM-41 silica pores disappears at pressure between 
1200 and 1600 bar. Xu et al. [38], by using simulations of TIP5P, ST2 water and other 
models for liquids with the LLCP, showed that the disappearing of the dynamic crossover 
of water can be associated to the presence of the LLCP (C). The simulations for water, 
indeed, display the same phenomenology of the experiments, with a crossover in the 
diffusion coefficient at P below the LL critical pressure P < Pc>, and with no crossover 
at P > Pqi. Xu et al. |38j presented numerical evidences that the crossover at P < Pc> 
is associated to the structural change occurring at the temperature of the maximum 
of the specific heat C™ ax along a line departing from the LLCP and extending in the 
one-phase region of the P-T phase diagram. This line, close to the LLCP, coincides 
with the Widom line [12], defined as the locus of the maximum correlation length in the 
one-phase region. 

Xu et al. [38] interpreted the absence of the crossover in the diffusion coefficient 
at P > Pqi as a consequence of the fact that the HDL-to-LDL spinodal occurs at 
almost constant P ~ Pqi for decreasing T. Therefore, when cooled at P > Pc, water 
never crosses the HDL-to-LDL spinodal and does not undergo the structural change 
responsible for the dynamic crossover. 

Hence, although the origin of the crossover has different interpretations [29], the 
experimental and numerical evidences suggest that the change in the dynamics is 
triggered by a local rearrangement of the H bond network |30j . Experiments and 
simulations, however, cannot give a definitive answer due to their finite resolution. 
For this reason is interesting to analyze the dynamic crossover of supercooled water 
within the framework of a Hamiltonian cell model that allows simulations and analytic 
calculations [TQl [12] . 

4. Water model 

We consider a cell model for water J9j [TUl [TTJ [12] based on the experimental observations 
that on decreasing P at constant T, or on decreasing T at constant P, (i) water displays 
an increasing local tetrahedrality [39l 140] . (ii) the volume per molecule increases at 
sufficiently low P or T, and (iii) the 0-0-0 angular correlation increases [13l [15] . 
consistent with simulations [22l I4T] . 

The system is divided into cells i & [1, . . . , N] on a regular lattice, each containing a 
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molecule, with a volume Vj > Vq, where Vq is the hard-core volume of one molecule, with 
a total volume V = Ylf v i- m ^ dimensions the distance between two nearest neighbor 
(n.n.) molecules i and j is r;j = [v\ + v l J d )/2. Since Vi is a continuous variable, the 
distance r it j is continuous. 

Every cell is occupied by a molecule. The dimensionless density for the molecule in 
cell i is v /vi G (0, 1]. We use a discrete two-state liquid- index ^ to quantify if the cell i 
is in the liquid phase or not, with rij = 1 if Vq/v > 0.5 and rii = otherwise. Therefore, 
J2f n i is the total number of liquid cells and iV — £\ n« is the total number of gas cells. 

The van der Waals attraction between the molecules is represented by the 
Hamiltonian term 

H = -e ^2 n i n i > ( 2 ) 
(id) 

where e > is the van der Waals attraction energy, which induces the liquid-gas phase 
transition. 

Each molecule % has four H-bond indices o~ij G [1, . . . , q], corresponding to four n.n. 
cells j, giving rise to q A different molecular orientations. Bonding and intramolecular 
(IM) interactions are accounted for by, respectively, the two Hamiltonian terms 

H B = -J^2 "i".i (S » ' ( 3 ) 

(id) 

where the sum is over n.n. cells, < J < e is the bond energy, 5 a ^ = 1 if a = b and 
&a,b — otherwise, and 

Him = -Ja n * S °ik°it> ( 4 ) 

i (k,£)i 

where Yl(ki)- denotes the sum over the IM bond indices (k,l) of the molecule i and 
J a > is the IM interaction energy with J a < J, which models the angular correlation 
between the bonds on the same molecule. 

When H bonds are formed, a small volume v-q is added to the local volumes Vi 
and Vj of the two H-bonded molecules i and j, increasing their average distance to 
r\j = (vi+Vj+VB) 1 ^ 1 2. Pictorially, this can be viewed as a local increase of the excluded 
volume associated with molecules i and j and is consistent with the experimental 
observation that H-bonded molecules form a low density open structure. Therefore, 
the total volume is proportional to the total number iVe of H bonds, as 

V = V + N B v B , (5) 

where Vq = Nvq is the volume of the liquid with no H bonds, and 

N B = Y n i n A^,^- ( 6 ) 

(id) 

Therefore, the cell model has the total Hamiltonian 

Htot = T~t + 7~Lb + 7~tiM (7) 
with the volume given by Eq. (jSJ). 
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5. Thermodynamics of the water model 

The cell model described in the previous section has been analyzed by mean field 
calculations M, HI H2l S21 S3] and Monte Carlo simulations [HI H21 S3] . The model's 
parameters considered in these studies are: J/e = 0.5, J a /e = 0.05, vb/v = 0.5 and 
q = 6. 

5.1. Mean field 

The mean field results [12], consistent with computer simulations [IT], reproduce the 
known phase diagram of fluid water, with the liquid-gas coexistence region ending in 
the critical point C, at ksTc/e = 1.03 ± 0.03 and PcVo/e = 0.18 ± 0.04, and with the 
temperatures of maximum density (TMD) at constant P that decreases with increasing 
P as in the experiments [HH [50] . 

In the deeply supercooled region the density has another discontinuity marking 
the coexistence region between two liquids at different densities. This discontinuity is 
associated with a discontinuity in m™ m , the value of the tetrahedral order parameter 
m a G [0, 1] that minimizes the molar Gibbs free energy g. Here, m a quantifies the 
number of H bonds with tetrahedral orientation and is defined as the difference between 
the number-density of o~ij = 1 and the average number-density for the other q — 1 states, 
i.e., m a = n a — (1 — n a )/(q — 1), where 



(iJ) 

and where the state = 1 corresponds to the appropriate state to form a (not 
bifurcated) H-bond with local tetrahedral order. 

In the NPT ensemble, the relevant free energy is the Gibbs energy per mole 



where u is the molar energy, s the molar entropy, and v the molar volume. As explained 
in Ref. [12], these quantities can all be written in the mean field approximation, giving 
rise to a mean field expression for g that can be minimized with respect to the gas-liquid 
order parameter m and the tetrahedral order parameter m CT . 

For any P at low k^T j e < 0.1, g has its minimum for m = 1 (liquid-phase value) and 
for a value m™ in that changes with T and P [Fig.flT])]. At constant P, m™ in decreases 
with increasing T, displaying a discontinuity above Pvo/e ~ 0.8 [Fig. (DO) right]. The 
discontinuity disappears at Pv^/e ~ 0.8 [Fig. ([I]) center] and the approach of m™ m to is 
always continuous below Pv Q /e = 0.8 [Fig.([T]) left]. The appearance of the discontinuity 
in m™ m denotes the occurrence at high P of a phase transition between two liquid 
phases with a different value of the tetrahedral order parameter m CT : an orientationally 
disordered phase {m a = 0) at high T, with no tetrahedral order, and a tetrahedrally 
ordered phase (m CT > 0.5) at low T. The phase separation disappears at the liquid- 
liquid critical point C, that in mean field is estimated at kBTc'/e = 0.062 ± 0.02 and 




(8) 



g = u — Ts + Pv , 



(9) 
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Tetrahedral Order Paramter m Tetrahedral Order Paramter m Tetrahedral Order Paramter m 

o a a 

Figure 1. The mean field molar Gibbs free energy g (dashed lines), in units of e, 
as function of the dimensionless tetrahedral order parameter m a 6 [0, 1] for different 
choices of T and P. The thick (red) line connects the points m™ ln of minimum g 
at different T for Pvo/e — 0.7 (left panel), 0.8 (center panel), 0.9 (right panel). In 
each panel the topmost line corresponds to the lowest ksT/e (0.06, 0.05 and 0.04, 
respectively) and the bottom line to the highest fc^T/e (0.08, 0.07, 0.06, respectively) 
with lines separated by ksST/e = 0.001. In all the panels m™ m increases when T 
decreases, being at the higher temperatures and ~ 0.9 at the lowest temperature. The 
value ra a = corresponds to the absence of tetrahedral order, i.e. to the high-density 
arrangement of water molecules. The value m a = 1 corresponds to full tetrahedral 
order of the H bonds, i.e. to the low-density arrangement of the water molecules. At 
Pv$/t = 0.7 (left panel), by increasing T, m™ m changes without discontinuity from 
~ 0.9 at ksT/e — 0.06 to at fcgT/e = 0.078, denoting a continuous change from 
the low-density arrangement at low T to the high-density arrangement at high T . At 
Pvo/e = 0.9 (right panel), instead, by increasing T, m™ m changes with discontinuity 
from ~ 0.6 at fcsT/e = 0.051 to at igT/e = 0.052, denoting a discontinuous 
phase change from the low-density liquid (LDL) to the high-density liquid (HDL). 
The pressure Pvo/e — 0.8 (center panel) is very close to the critical pressure Pc, 
below which the discontinuity seen at higher P disappears. 

Pc v o/ e — 0.82 ± 0.02 [12]. The discontinuity in m a determines the discontinuity in the 
density between the two phases, as can be shown analytically p2], separating the liquid 
in LDL at low T and HDL at high T, consistent with the LLCP scenario. 

5.2. Monte Carlo 

In the Monte Carlo simulations of the cell model, the Hamiltonian Titot in Eq.Q is 
simulated with the standard spin-flip method [TTl |4*2| |4"3] . For sake of simplicity, we 
solve the model under the condition of being in a homogeneous phase with v i = v for 
any i. Although this condition can be easily removed, the solution of the model under 
this condition gives good qualitative agreement with the experiments, as we will see in 
the following sections. 

The model is defined in any dimension. However, since we assume for simplicity 
that any molecule can form at most four H bonds, we solve the model on a regular 
square lattice of cells. This choice is particularly appropriate for the study of quasi-2D 
water [44] between hydrophobic surfaces and is a first order approximation to the layer 
of water between the surface of a protein and the surface of the frozen bulk water at low 
T. Of course, in this approximation the interaction with the confined surfaces is not 
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taken into account [45J. Nevertheless, the qualitative comparison with the experiments 
is satisfactory. The lack of the third dimension, and the number of neighbors limited 
to four, could affect properties related to the amount of free space around a molecule, 
such as the diffusion constant [HH H7J HH] . These properties are not considered in the 
analysis of the model reported here. 

To further simplify the simulations, the van der Waals interaction, represented in 
mean field by the Hamiltonian term in Eq.(j2J), is replaced by an equivalent Lennard- 
Jones potential with characteristic energy — e at distance Rq = v\ and truncated with 
a hard-core in its minimum at Rq [11] H2J H3]. Simulations performed with a standard 
Lennard- Jones potential, without truncation, show that the results are qualitatively 
unchanged |47J. 

The resulting phase diagram [TT] displays, at high T and low P, a gas-liquid phase 
transition ending in a gas-liquid critical point C. From C departs in the supercritical 
region the gas-liquid Widom line, i.e. the line of maximum — but finite — correlation 
length for the fluid. At lower T, the phase diagram displays the TMD line, retracing 
toward low T at high P as in the experiments [3]. At lower T and high P, the phase 
diagram displays a HDL-LDL phase transition with negative slope in the P-T plane, 
ending in a LL critical point C. As for the gas-liquid critical point, also from C the 
line of maximum correlation length for the liquid, the LL Widom line [38j [12] , departs 
into the one-phase region (Fig. EJ). Therefore, the Monte Carlo results confirm the mean 
field analysis, displaying a phase diagram consistent with the LLCP scenario. 

Recent calculations p2] show that by varying the parameters of the cell model, it is 
possible to recover also (i) the singularity free scenario, (ii) the stability limit scenario 
and (Hi) the critical point free scenario. In the following we report about the analysis of 
the dynamics of the model for the LLCP scenario and the singularity free (SF) scenario. 



6. Dynamics of the water model 

We first consider the case of the LLCP scenario. Performing Monte Carlo simulations 
for T and P around the values of the LL Widom line T\y(P), Kumar et al. [12] found 
that the correlation time r of Si = J2j vl^i which quantifies the degree of total bond 
ordering for site i, displays a dynamic crossover from Vogel-Fulcher-Tamman (VFT) 
function at high T to Arrhenius T-dependence at low T (Fig. [3ja). The VFT function 
is given by 

T-T . 

where r^ FT , T\, and T are three fitting parameters, and the Arrhenius function is given 
in Eq. (prj). This result is qualitatively consistent with the dynamic crossover found in 
experiments for confined water and hydration water [231 [2E1 [32], [33], [36] [37] and has been 
related to the presence of the LLCP [Ml 154] . 

It is, therefore, interesting to compare this result with the dynamic behavior of 
the cell model in the hypothesis in which the SF scenario holds. To do this, Kumar et 



VFT _ VFT 

t = r exp 



(10) 
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Figure 2. The Monte Carlo phase diagram for the cell model at low T for N = 3600 
water molecules [13] ■ C is the HDL-LDL critical point, end of first-order phase 
transition line (thick line) [5]; symbols are maxima of the coefficient of thermal 
expansion ap| max (O)j isobaric specific heat C™ ax (□), \dps / 'dT| max (O) the numerical 
derivative of the probability of forming a H-bond, proportional to the fluctuation of 
the number of bonds (<5 2 A^B) max (A); the Widom (solid) line, corresponding to locus of 
maximum correlation length and estimated as the average between |ap| max and Cp ax , 
coincides within the error bars with \dp^/dT\ max 1 i.e. with the locus of the maximum 
structural variation. Dashed lines are guide for the eyes. 



al. [42J analyzed the dynamics of the cell model when the H-bond correlation is zero 
( J rj = 0) and the cell model recovers the SF model of Sastry et al. [20] . 

The result (Fig. Ob) [42] shows that also in this scenario a dynamic crossover is 
expected. The temperature of the crossover is coinciding, within the numerical precision, 
with the temperature T(C™ ax ) of maximum isobaric specific heat, which in the SF 
scenario play a role equivalent to 7\y(P) of the LLCP scenario. Indeed, both T(Cp ax ) 
and T W (P) mark the temperature T^ruct °f the maximum structural change for the 
liquid [43] . 

The cell model allows to clarify that the dynamic crossover is, indeed, a direct 
consequence of the structural change occurring at T™* ct . By calculating, in the mean 
field approximation, the T-dependent activation energy Ea{T) mf necessary to (z) break 
a non-tetrahedral H bond, (ii) reorient the water molecule and (Hi) form a new H bond 
in a tetrahedral orientation, Kumar et al. [43J calculated the mean field correlation time 
as 

'E A (T) MF ~ 



r MF = r exp 



, T . (11) 
k B T 

and compared r MF with the correlation time from the Monte Carlo simulations, finding 
an excellent agreement. Since Ej^ F is an explicit function of the number Nb of H 
bonds, therefore r MF is a function of the structure of the liquid |43j. Since iVe rapidly 
changes for T > T™^, has its maximum variation at T™* ct , and slowly changes for 
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Figure 3. Monte Carlo results for the dynamic crossover in the orientational relaxation 
time r of the cell water model for a range of different pressures [42 . (a) In the case of 
the liquid- liquid critical point (LLCP) scenario , the crossover occurs at a temperature 
consistent with the Widom line 7V(P). (b) In the case of the singularity free (SF) 
scenario, the crossover occurs at a temperature consistent with the T(C™ ax ). In both 
panels the large hatched circles mark Tw{P) and T(Cp ax ), with a radius approximately 
equal to the error bar. Solid and dashed lines represent Arrhenius and Vogel-Fulcher- 
Tamman fits, respectively. In both cases, Kumar et at [42j predicted that the dynamic 
crossover occurs at approximately the same value of r for all the values of pressures 
studied, i.e. the crossover is isochronic. 



T < T™* ct , the correlation time is a no n- Arrhenius function of T for T > T™* ct and 
almost an Arrhenius function for T < T s ™^ t , giving rise to a dynamic crossover [43J. 
The agreement between the mean field calculation and the Monte Carlo results clarifies 
that the low-T regime is only approximately Arrhenius. 

The analysis of the cell model allows, furthermore, to make a number of predictions 
about the dynamics of water at low T. 

(i) At the dynamic crossover the value of r is about the same for all the pressures 
considered in Ref. [42], i.e. the crossover is isochronic and the correlation time at 
the crossover is t~l = 10 3 / 2 in units of Monte Carlo steps [12] ■ Comparison with 
experiments, in the next section, allows to convert this result in real time units. 

(ii) The activation energy Ex, that is found by fitting the low-T dynamic regime with 
an Arrhenius behavior, decreases linearly for increasing P [42] . 

(iii) If we fix a characteristic time, for example the crossover time tl, the temperature 
Tl at which this correlation time is reached decreases linearly for increasing P |42j . 

All these predictions are verified in both the LLCP scenario and the SF scenario. 
However, Kumar et al. [12] found a difference between the two scenarios: for the LLCP 
scenario the index Ex/ '(/cb^l) increases for increasing P, while for the SF scenario 
Ex/iksT-if) is constant. However, the predicted increase of -Ea/(&b2l) is of the order 
of 1% [12] ■ As we will discuss in the next section, this increase is within the present 
experimental error bar. 
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7. Comparison with the experiments on hydration water 

Recently, Chen and coworkers have performed extensive QENS experiments to 
investigate the dynamical behavior of hydration water on lysozyme, spanning a range 
of pressures going from ambient pressure up to 1600 bar [5T] . By measuring the mean 
square displacement of the of H-atoms in lysozyme, they found a dynamic crossover for 
all pressures studied [51] . They also found that the translational correlation time of the 
H-atoms of the hydration water shows a dynamic crossover for P — 1, 400, 800, 1200, 
1500 bar [51]. For P < 1500 bar, they showed that the lysozyme dynamic crossover 
and the hydration water dynamic crossover occurs at the same T [51], extending their 
previous results on lysozyme [23], DNA [32] and RNA [33] at ambient pressure. 

They also showed that the dynamic crossover for the hydration water disappears 
at P = 1600 bar [51J. This finding, together with their previous measurements for 
water confined in MCM-41 Silica pore [37], suggests that the disappearing of the non- 
Arrhenius to Arrhenius crossover at P = 1600 bar is independent of the constraint used 
to avoid homogeneous ice nucleation. As shown in Ref. [38], the disappearing of the 
crossover could be the consequence of the LLCP occurring at a pressure between 1500 
and 1600 bar. 

The comparison of the experimental data in Ref. [51] with the results of the Monte 
Carlo simulations allows several observations. 

• The correlation time tl at the crossover is constant as predicted by the cell model, 
point (i) in the previous section. 

• r L = 10 3//2 in units of Monte Carlo steps (MCS) [32] corresponds to 10 4,25 ps, i.e. 
1 MCS^ 600 ps (Fig. H. 

• The high-T (non-Arrhenius) behavior of the correlation time in real water has a 
stronger dependence on T than the one seen in the cell model at P = (Fig. HJ) and 
is more similar to the behavior seen in the model at higher P. Since in the model 
this behavior is regulated by the increase of number N-g of H bonds for decreasing 
T, this observation implies that in real water rapidly increases approaching the 
temperature of maximum structural change T™ r ^* t , and is smaller in real water 
than in the model for T > T^ x ct [33]. 

• The low-T (approximately Arrhenius) behavior is well reproduced by the cell model 
(Fig. HI). This result implies that the estimate of the activation energy at P = in 
the model compares with a reasonable agreement to the activation energy measured 
for hydration water. 

• The direct comparison of the Monte Carlo estimated activation energy Ea and the 
QENS data (Fig. [5J upper panels) shows that the prediction (ii) of the previous 
section is verified in the experiments, with Ex decreasing linearly for increasing P. 
The agreement is only qualitative, since the variation of E\ in the cell model is too 
large with respect to the experiments. 



Pressure Effects in Supercooled Water: Comparison between a 2D Model of Water and Experiments for Surf 



P = 1bar 
P = 400 bar 
P = 800 bar 
P = 1200 bar 
(MC) P = 




0.002 



0.004 0.006 
1/T(K" 1 ) 



0.008 



Figure 4. Comparison between the Monte Carlo [H] results and the QENS data [ST] 
for the correlation time r of supercooled water. Monte Carlo results (open symbols) 
are for P = and are rescaled in such a way that the crossover occurs at at the same T 
and r of the experimental data (full symbols). The experimental data are for pressures 
going from ambient to 1200 bar. 



• The crossover temperature /c_bTl in the experiments verifies the prediction (iii) of 
the cell model, i.e. k B T L decreases linearly for increasing P (Fig. [SJ lower panels). 
Also in this case, the variation in the model is too large and the agreement is only 
qualitative. 

• The crossover temperature /cbTl/c ~ 0.16 at ambient pressure corresponds to 220 K 
(Fig. HJ), leading to a van der Waals interaction strength e ~ 11 kJ mol -1 and 
a H-bond strength J = 0.5e = 5.5 kJ mol -1 . This H-bond strength is about 
1/4 of the H-bond strength expected below 100° C [52J and is about the value 
estimated for then van der Waals attraction based on isoelectronic molecules at 
optimal separation [53J. 

We finally observe that the prediction of the model about the different behavior of 
the index E\/ (ZcbTtJ in the two scenarios, cannot be verified in the experiments, since 
the predicted difference between the two scenarios (of the order of 1%) is within the 
error bars of the measurements [51] . The QENS data in Ref . [51] show a non-monotonic 
behavior for this index within an error bar larger than the 1% variation predicted for 
the LLCP scenario. Hence, the data could be considered also constant within the error 
bars, as predicted for the SF scenario. Therefore, is not possible to discriminate between 
the two scenarios on the basis of the present data for the index Ea/ {ksT^). 
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Figure 5. Monte Carlo results and the QENS data [5T] for the activation energy 
Ea of the low-T regime and the crossover temperature Tl, as functions of the pressure 
P. Upper left panel: Monte Carlo results for Ea [32] : for sake of comparison with 
the experimental data, pressure is rescaled by an arbitrary factor 1500 bar/(0.8e/uo) 
and energy by a factor 16.4 kcal mol _1 /e ~ 68.7 kJ mol~ 1 /e; both the liquid-liquid 
critical point (LLCP) scenario results (circles) and the singularity free (SF) scenario 
results (squares) are reported. Lower left panel: Monte Carlo results for the crossover 
temperature TL; pressure is rescaled as above and temperature is rescaled by a factor 
220 K/(0.16e/fcs); symbols are as in previous panel. Upper right panel and lower right 
panel: QENS data for the lysozyme hydration water [51]. 

8. Discussion and Conclusions 

Several questions remain open. The work reported here allows us to propose possible 
answers to some of these questions. 

• Why does water have anomalies? 

Simulations [6] first and, after, experiments [13] have shown the existence of two 
local configurations for water molecules: tetrahedral up to the second shell (or open) 
and non-tetrahedral (or closed). These local arrangements are the consequence of 
the orientational character of the H bonds and their typical competition between 
attraction and repulsion. Many of the anomalies can be understood on the basis 
of the tendency of the H bond to form open configurations to minimize the energy, 
frustrated by the necessity of reducing the occupied volume at high density or 
pressure. 
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• What are the implications of these open and closed configurations? 

The cell model discussed here [TO], [TTJ, [12] shows that the existence of these 
local configurations, related to the tendency of the H bonds to correlate and order 
in a tetrahedral way, is enough to imply the occurrence of a liquid-liquid phase 
transition, possibly ending in a critical point. Only in the hypothesis that the H 
bonds formed by the same molecule are completely uncorrelated, as in the Sastry et 
al. model [20) . the open and closed configurations do not separate into two phases, 
giving rise to the singularity free scenario. 

• Why is there a dynamic crossover in the correlation time of supercooled water? 
Experiments [23J and simulations [38] offer evidences that the dynamic crossover 
is due to the local variation of H bond network [34|. The cell model discussed 
here, furthermore, offers an analytic relation, based on a mean field approximation, 
between the dynamic crossover and the structural change [32J [12]. In the cell 
model it is assumed that the relevant mechanism for the low-T water dynamics is 
the breaking of the H bonds to reorient the molecule and form more tetrahedral 
bonds, as seen in simulations [33] . 

• Why is the dynamic crossover at the locus T(C™ ax ) of maximum specific heat? 

At T(Cp ax ), the number N-q of H bonds has its maximum variation. The liquid 
changes from the HDL-like local structure at high T to LDL-like local structure at 
low T. In the HDL-like regions Nb rapidly increases when T decreases, implying 
that the activation energy Ex for the breaking and reorientation of the H bonds 
increases and the behavior is non-Arrhenius [43J. In the LDL-like regions, 
is almost constant with T, implying that Ex is approximately constant and the 
behavior is asymptotically Arrhenius [43J. 

• Is the dynamic crossover an evidence of the liquid-liquid critical point? 

The dynamic crossover is consistent with both the liquid-liquid critical point 
scenario and the singularity free scenario [32]. However, in the hypothesis of the 
presence of a liquid-liquid critical point is possible to rationalize the disappearing 
of the dynamic crossover at high pressure [371 [3H [51], while in the other scenarios 
no change in the crossover is expected with the increase of P. 

• Is the dynamic crossover a cooperative process? 

The dynamic crossover is related to the increase of the correlation length £ in the 
liquid. This quantity increases for increasing P along the Widom line and diverges 
at the liquid-liquid critical point. Therefore, the dynamic crossover is a cooperative 
process at the liquid-liquid critical point. 

• How pressure studies can help us in understanding the physics of water? 

The recent experiments on protein hydration water [51] show that is possible to 
analyze the dynamic crossover of water at very low T, measuring quantities, such 
as the crossover temperature T L and the activation energy Ex, whose behaviors 
have been predicted by the theory [42] . The experiments have verified three of the 
four theoretical predictions, being the forth within the error bars. 
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In conclusion, in this report we compare the predictions of a cell model for water 
analyzed by MC in 2D [42] with the QENS data for water at the surface of lysozyme 
[51~] . The data show a dynamic crossover from high-T non-Arrhenius behavior to low-T 
quasi-Arrhenius behavior for the water relaxation time at P < 1600 bar. Both the 
temperature 7l(P) and the low-T average activation energy E\(P) linearly decrease for 
increasing P, as predicted by the model. The relaxation time at the crossover t(Tl) is 
independent of P, as in the model. It is an open question if this isochronic behavior is 
related to the constant-time dynamic crossover observed in glass formers [55] . 

The mean field approach [42, 43j allows to find a functional relation between E\ and 
the number Nb(T,P) of H bonds, i.e. allows to show in a clear way that the dynamic 
crossover is a consequence of the structural changes occurring at the temperature of 
maximum specific heat T(Cp ax ). However, the comparison with the experiments does 
not clarify which scenario describes better the supercooled region of water. Indeed, 
within the experimental error bar larger than 1%, both the SF and the LLCP scenario 
could be consistent. Nevertheless, the disappearing of the crossover observed in the 
experiments at P = 1600 bar is expected only in the presence of a LLCP. 
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